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^ ■ Abstract A gravitational wave must be nonlinear to be able to transport its own source, 

£NJ , that is, energy and momentum. A physical gravitational wave, therefore, cannot be represented 

O \ by a solution to a linear wave equation. Relying on this property, the second-order solution 

describing such physical waves is obtained. The effects they produce on free particles are found 
to consist of nonlinear oscillations along the direction of propagation. 
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1 Introduction 

For a long period after the advent of general relativity, the question of the existence or 
^jjy not of gravitational waves was a very controversial issue. In the seventies, the discovery 

of a binary pulsar system whose orbital period changes according to the predicted wave 
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emission put an end to the controversy pQ. In fact, that discovery provided an indirect 
but compelling experimental evidence for the existence of gravitational waves [3] . That 
evidence, however, did not provide any clue on their form and effects. As a matter of fact, 
although widely considered a finished topic [3] , it actually remains plagued by many obscure 
£Sj ■ points [5j- For example, although there seems to be an agreement that the transport of 

, energy-momentum by gravitational waves is essentially a nonlinear phenomenon, instead of 

going to the second order, one usually assumes a "mixed" procedure, which consists basically 
in assuming that gravitational waves carry energy (are nonlinear, or at least second order), 
but at the same time, because this energy is very small, one also assumes that its evolution 
can approximately be described by a linear (first order) equation [5J. When one speaks 
of "linear gravitational waves", therefore, one means nonlinear gravitational waves whose 
dynamics is assumed to be approximately described by a linearized equation. This means 
that, in addition to the sequential levels of accuracy implied by the perturbative analysis, 
there is also another approximation, implied by the "mixed" approach, according to which 
all first-order equations describing a gravitational wave are to be interpreted as only nearly 
correct [7]. 

Such assumption, however, is an unjustified surmise: the issue is not a matter of ap- 
proximation, but a conceptual question A gravitational wave either does or does not 
carry energy: if it carries, no matter how small it is, it cannot satisfy a linear equation. 
It is, therefore, conceptually unsatisfactory to assume that a gravitational wave satisfying 
a linear equation is able to transport energy and momentum. If applied to a Yang-Mills 
propagating field [9], the approximation described above would correspond to assume that, 



1 There are other arguments against this assumption. See, for example, Ref. [8]. 
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for a field with small-enough amplitude, its evolution can be accurately described by a linear 
equation. Of course, this is plainly wrong: a Yang-Mills propagating field must necessarily 
be nonlinear to carry its own (color) source, otherwise it is not a Yang-Mills field. Analo- 
gously, a gravitational wave must necessarily be nonlinear to transport its own source — that 
is, energy and momentum. 

Taking into account these premises, a critical review of the gravitational wave theory has 
been published recently [TU] . In that paper, it was discussed why the standard approach to 
the gravitational wave theory is not satisfactory. Here, instead of using the mixed approach, 
we proceed to the second order and obtain the corresponding nonlinear gravitational wave. 
It is important to remark that this re-interpretation of the gravitational wave concept has 
no implications for the usual expressions of the power emitted by a mechanical source. 
In particular, the (nonlinear) quadrupole radiation formula gives a correct account of the 
energy emitted by a binary pulsar, for example. The only we claim is that the energy and 
momentum are not transported away by linear, but by nonlinear waves. The basic purpose 
of the present paper is to make an analysis of these waves, as well as of their effects on test 
particles. 



2 Linear Approximation 

2.1 Linear Wave Equation 

The study of gravitational waves involves basically the weak field approximation of Einstein 
equation 

where Q^ u is the source energy-momentum tensor. That is arrived at by expanding the 
metric tensor according to 

V^cT = rT +eh% + e 2 h% + ..., (2) 

where e is a small parameter introduced to label the successive orders in this perturbation 
scheme. When the metric tensor is expanded according to ([2]), we are automatically as- 
suming that there is a background Minkowskian structure in spacetime, with metric rj^ u . 
Accordingly, the gravitational waves are interpreted as perturbations 

h» u = eh% + £ 2 h% + ... (3) 

propagating on that fixed Minkowskian background. This interpretation is consistent with 
general relativity, as well as with the point of view of field theory, according to which a field 
always propagates on a background spacetime 

Assuming expansion ([2]), the first order Ricci tensor is 

R(T)ij,u = ^xY^fjj, — d u T^ r ^x . (4) 
Using the first order Christoffel connection 

r^j/w = \ (^~d^h^ u — dvh^n + d^h(x)jj,v + d^h^ + d u h^ — \f\ ilv d^h (1 ^j , (5) 
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the Ricci tensor and the scalar curvature become, respectively, 

R(Dtiu = 5 ( n /i(i )M v - lrifj, v ah m - dxdph^v - d u d x h (1)pX ^j (6) 

and 

R w = -±nh {1) -d x d p h%, (7) 

where □ = rj px d p d x is the flat spacetime d'Alembertian, and h w = h^x- in consequence, 
the first order sourceless gravitational field equation becomes 

a h Wfll/ - dxdfji^v - d u d x h {1)pX + V^dxdph^ = 0. (8) 

Now, as is well known, wave equation ([8]) is invariant under (infinitesimal) general 
spacetime coordinate transformations. Analogously to the electromagnetic wave equation, 
which is invariant under gauge transformations, the ambiguity of the gravitational wave 
equation can be removed by choosing a particular class of coordinate systems — or gauge, 
as it is usually called. The most convenient choice is the class of harmonic coordinate 
systems, which at first order is fixed by 

dphft = 0. (9) 
In this case, the field equation ([8]) reduces to the relativistic wave equation 

□ h p (1)P = 0. (10) 



2.2 Linear Waves 

A monochromatic plane-wave solution to equation (|l(jp has the form 

fyi)AtK = A {1)flu exp[ik p x p ], (11) 

where A Wpu = Am™ is the polarization tensor, and the wave vector k p satisfies 

k p k p = 0. (12) 

The harmonic coordinate condition ([9]), on the other hand, implies 

k P hf 1)v = 0. (13) 

In analogy with the Lorentz gauge in electromagnetism, it is possible to further specialize 
the harmonic class of coordinates to a particular coordinate system. Once this is done, the 
coordinate system becomes completely specified, and the components A Wpu turn out to 
represent only physical degrees of freedom. A quite convenient choice is the so called 
transverse-traceless coordinate system (or gauge), in which [12] 

h p {1)P = and h^ v U^=0, (14) 

with UY 0) an arbitrary, constant four-velocity. 
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Now, although the coordinate system {x^} has already been completely specified (the 
transverse-traceless coordinate system), we still have the freedom to choose different local 
Lorentz frames e a . Since the metric g^ u = n a } ) e a ^e h v is invariant under changes of frames, 
the metric perturbation will also be invariant. In particular, it is always possible to choose 
a specific frame, called proper frame, in which UY^ = 5q u . In this frame, as can be seen from 
the second of the Eqs. ([ED , 

fc&o = (15) 
for all fx. Linear waves satisfying these conditions are usually assumed to represent a plane 
gravitational wave in the transverse-traceless gauge, propagating in the vacuum with the 
speed of light. Its physical significance, however, can only be determined by analyzing the 
energy and momentum it transports. 

2.3 Energy and Momentum Transported by Linear Waves 

The energy-momentum tensor of any matter (or source) field tp is proportional to the func- 
tional derivative of the corresponding Lagrangian with respect to the spacetime metric. 
Since such a derivative does not change the order of the Lagrangian in the matter field tp, 
both the Lagrangian and the energy-momentum tensor will be of the same order in the 
field variable tp. For example, both Maxwell's Lagrangian and its corresponding energy- 
momentum tensor are quadratic in the electromagnetic field. Now, it is a well known fact 
that the gravitational field is itself a source of gravitation. This means that the gravita- 
tional energy-momentum current should appear explicitly in the gravitational field equation. 
Accordingly, the wave equation (jlOp should read 

\QttG 

D h WlJLU = ^ 4 t (l)liu . (16) 

At the linear approximation, therefore, the gravitational energy-momentum density t m ^ u is 
restricted to be linear. However, since the energy-momentum density is at least quadratic 
in the field variable, t WfJjl/ vanishes in the linear approximation, leading to the wave equa- 
tion (fTUjl . 

The above property is a crucial difference between linear gravity and electromagnetism, 
and is often a source of confusion. Even though the electromagnetic waves are linear, they 
do transport energy and momentum. There is no any inconsistency in this result because 
neither energy nor momentum are sources of electromagnetic field, and consequently the 
energy-momentum tensor does not appear explicitly in the electromagnetic field equation. 
In other words, even though the electromagnetic field equations are linear, the energy- 
momentum tensor is not restricted to be linear. The linearity of the electromagnetic wave 
equation, however, restricts the electromagnetic self-current to be linear, and consequently 
to vanish. This means that the electromagnetic wave is unable to transport its own source, 
that is, electric charge. A linear gravitational wave is similarly unable to transport its own 
source, that is, energy and momentum. Only a nonlinear wave will be able to do it. This a 
subtle, but fundamental difference between electromagnetic and gravitational waves. 

The consistency of this result can be verified by analyzing the generation of linear waves. 
In the presence of a source, the first order field equation reads 

16vrG 

Q h(l)flV — -J (17) 
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with @(i)uu the first order source energy-momentum tensor. As a consequence of the coor- 
dinate condition Q, it is easy to see that 



d"Q Wl „, = 0. (18) 

Instead of the usual covariant derivative, is conserved with an ordinary derivative at 

the first order. Since this is a true conservation law, in the sense that it leads to a time 
conserved charge, we can conclude that in the linear approximation a mechanical system 
cannot lose energy in the form of gravitational waves o As discussed in section 1, this 
problem is usually circumvented by assuming the mixed approach, according to which this 
equation is to be interpreted as nearly true. 



2.4 Generation of Linear Waves 

Let us consider the first order field equation (|17p . A solution is the retarded potential 

_ AG f d 3 x' , 

with the source considered in the retarded time 

t> = t _ W-x'\ ^ (2Q) 



c 

At large distances from the source we can expand 

\x — x'\ ~ r — x' ■ n + . . . , (21) 

where r = \x\ is the distance from the source, and h is a unit vector in the direction of x. 
The leading order term of hmuv is obtained by replacing \x — x'\ in the denominator of 
Eq. ([19]) with r, 

h(i)fj,u = ^4 J d?x @ (1 ) /ii y(t , x ), (22) 

where now 

. r x' ■ h , . 

t' = t-- + . (23) 



The Fourier transform of 8(i)m/ is 

©cWf , ^) = / ®vA", k) e-^-Hft*. (24) 

Substituting in Eq. (|22p and performing the integrations in d?x' and d 3 fc, we obtain [2] 

hwv» = — i r d ^ e (1)F (^n/c)e-^). (25) 



o 



2 This is consistent with the fact that linear gravitational waves do not transport energy nor momentum. 
The existence of a linear solution is a mere consequence of the use of a perturbative scheme, but alone it 
does not represent the physical wave. 
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We see from this expression that, if the source oscillates with a single frequency oj, the 
plane wave h Wpi/ will necessarily propagate with the same frequency. However, we know 
from the quadrupole radiation formula that, if the source oscillates with frequency oj, the 
gravitational radiation should come out with frequency 2wH The reason for this factor of 
2 is that both the generation and the effects of gravitational waves on free particles are 
essentially tidal effects, which we know to occur twice during a complete cycle. This is a 
clear indication that h^p U alone cannot represent the physical gravitational wave. 



3 Second Order Approximation 

3.1 Second— Order Wave Equation 

At the second order of the iterated perturbation scheme, the harmonic coordinate condition 
reads [TT] 

dphhv = 0. (26) 

In these coordinates, the second order gravitational field equation can be written in the 
form 

_ 167rg ( t P , p \ (27 \ 

<-> "(2)f — ^4 V L {2)V T ^J(1)V) J \^'J 

where t p 2)U = t p 2)V {h w , /i {1) ) represents all terms coming from the left-hand side of Einstein 
equation, in addition to the d'Alembertian term. It can be interpreted as the second order 
energy- momentum pseudotensor of the gravitational field [6j. 

Far away from the sources, the second order gravitational waves are governed by the 
sourceless version of the wave equation (l2~T1). 



nh^ = ^-l^ = N^(h ll) ,h il) ), (28) 
where, already considering the traceless gauge condition h w = 0, 

N^(h w ,h w ) = -h%d p d h% + ±d»h Wpi ,ff'h% + d a h^{dPh v {l)p + d p h%) 

-dfhvprdPhfi - PhvpadPhft + ^ (dph^rd^h^ - \d T h Wpa d T h^ . (29) 

Using for hrX the plane wave solution (jlip . the wave equation becomes 



h u A pa — -k^h" A A pa — k k a A pp A v 
h h A^P A va 4- h u hP A AP a 4- hP-hP A A va 



if ( r 1 

2~y^p^ ^wcrT^fi) — —k T k A^p a A^ 



Use of the constraints (fT2l) and (PT31) reduces it to 



) 

exp[i2kpx p }. (30) 



3 See, for example, Ref. [2], page 105. 
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with 

* (9) = A {1)pa A%. (32) 

It is worth mentioning that the second-order wave equation is quadratic in the first-order 
solution h^ v . The factor "2" in the exponential of the right-hand side is a reminder of this 
nonlinear, quadratic dependence. 



3.2 Second— Order Nonlinear Waves 



A general solution to the wave equation (|3ip is given by a solution to the homogeneous 
equation plus a particular solution to the non-homogeneous equation. A monochromatic 
traveling-wave solution can then be written in the form 



K2) = i A (2) + iB m ) ex P^ 2k p xP ]> 



A {2) 



16 ' 



(33) 
(34) 

8 K a k a " ' ^ 

with K a an arbitrary wave number four-vector. As a direct inspection shows, this solu- 
tion satisfies the harmonic coordinate condition (|26|) . The physical gravitational wave is 
represented by the real part of the solution, that is, 



where 



and 



n (2) 



$(2) Kgx b 



k^k v 



h% = A% cos[2k p x p ] - B% sm[2k, 



x'- 



(36) 



Observe that the amplitude BtfZ depends explicitly on the wave number — or equivalently, 
on the frequency of the wave. This is a typical property of nonlinear waves. 
The amplitude of the first part of the solution satisfies 



AH 



A 



(2) 



(2) 



1 



and k^A^ — —k v A^ 2 ). 



(37) 



We consider now a laboratory frame — with a Cartesian coordinate system — from which 
the wave will be observed. In this case, only the diagonal components of A™ are non- 
vanishing and obey 



More specifically, 



A a 



A yy = A z 



A . 



(38) 





A 
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The second part, on the other hand, satisfies 

B 



- D (2)M 



'(2) 



and kfj,Bj% 



0. 



(39) 



(40) 
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If we consider, for example, a wave traveling in the z direction of the Cartesian system, for 
which 

k p = (w/c, 0, 0, o;/c) , (41) 
the coefficient will be of the form 



8K n k a c 2 











1\ 


























M 








V 



(42) 



Considering both parts of the solution we see that the second order wave is neither transverse 
nor traceless. 

As already seen, if r denotes the distance from the source, the amplitude of the first 
order solution scales according to A^K ~ 1/r. As an immediate consequence, $( 2 ) ~ 1/r 2 . 
This means that the amplitude of the first part of the solution (|36p falls off as 

~ 1/r 2 . (43) 

Due to an additional linear dependence on the distance, the amplitude of the second part 
falls off as 

B% ~ 1/r. (44) 
At large distances from the source, therefore, the dominant solution will be of the form 

^~£^sin[2Va (45) 

Usually, second-order effects are supposed to fall off as 1/r 2 , and for this reason they are 
assumed to be neglectful at large distances from the source [12]. However, as shown above, 
the second-order gravitational wave h^L falls off as 1/r, and consequently the arguments 
used to neglect them are not valid in this case. Observe also that, if the source oscillates 
with a single frequency u, the field will propagate, as appropriate for a quadrupole 
radiation, with a frequency 2ui [2]. This factor of 2 is a direct consequence of the nonlinear 
nature of the gravitational wave (see the comment just below Eq. (|32p ). and provides one 
more evidence that — and not — represents the physical (quadrupole) gravitational 
wave. 

3.3 Generation of Nonlinear Waves 

As can be seen from Eqs. (|26|) and (|27|) . the second order total energy-momentum tensor is 
conserved: 

+ %A = °- ( 46 ) 

The source energy- momentum tensor, on the other hand, as determined by the second order 
Bianchi identity, is conserved only in the covariant sense: 

v (l"(2)i' — u p KJ (2)V T -L (i)p/j, ^(1)^ 1 (1)V[± w (l)p — u - I 

At the second order, therefore, the source energy-momentum tensor is not truly conserved — 
it does not lead to a conserved charge. As a matter of fact, the above covariant conservation 
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law is not a true conservation law, but simply an identity governing the exchange of energy 
and momentum between gravitation and matter [15]. As a consequence, in contrast to what 
happens at the first order, at the second order a mechanical system can lose energy in the 
form of gravitational waves. 

It is important to remark once more that the usual expressions of the power emitted 
by a mechanical source, and in particular the quadrupole radiation formula, give a correct 
account of the energy emitted by a mechanical system. The reason is that nonlinear methods 
have always been used in the study of wave generation by such systems. Furthermore, 
the quadratic energy- momentum pseudotensor is the complex traditionally used to 
calculate the energy and momentum transported by gravitational waves. What we claim 
here is that, instead of being transported by the linear waves h^X, this energy is actually 
transported by the second-order gravitational wave h^- Notice from Eq. (|27p that 
appears as source of the second-order gravitational field h^ 2 y It represents, therefore, the 
energy and momentum transported by the second-order gravitational waves. 



4 Effects on Free Particles 

4.1 The Geodesic Deviation Equation 

Let us consider, as usual, two nearby particles separated by the four- vector £ a . This vector 
obeys the geodesic deviation equation 

Vf/Vt/r = R a ilv pU»U v f, (48) 

where = dx^/ds, with ds = g^ u dx^dx u , is the four- velocity of the particles. Now, each 
order of the gravitational field expansion 

R a fiu/3 = £ + £ 2 R%)iivp + • ■ • , (49) 

which follows naturally from ([2]), will give rise to a different contribution to £ a . For consis- 
tence reasons, therefore, this vector must also be expanded: 

r = Co)+^u>+^(% + --- • ( 5 °) 

In this expansion, represents the initial, that is, undisturbed separation between the 
particles. As the four-velocity depends on the gravitational field, it should also be 
expanded. However, since the gravitational wave is interpreted as a perturbation of the flat 
Minkowski spacetime, the movement produced on free particles will also be considered in 
Minkowski spacetime. This means that we can write = Ufa-, = dx^/ds^, where 

dsfo) = V^u dx^dx u (51) 

is the flat spacetime quadratic interval. Of course, the four-velocity depends on the 
choice of the initial condition — or equivalently, on the choice of the local Lorentz frame 
from which the phenomenon will be observed and measured. Using then the freedom to 
choose this frame (see section I2.2() , we can choose a frame fixed at one of the particles - 
called proper frame. In that frame, the proper time s (0) coincides with the coordinate x° [7], 
and the particle four-velocity assumes the form 

U{ 0) = 6» = (1,0,0,0). (52) 
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4.2 First-Order Effects 

Considering that £^ represents simply the undisturbed separation between the particles, 
at the lowest order the geodesic deviation equation is 



ds 2 



+ Ufa) d p (r" )( g 7 U( 0) ) £( 0) — R^pi/p ?7(o) ^(o) £(o)- (53) 



'(0) 

Substituting £7/1 as given by Eq. ([52]) . we get 



ds 2 



- + <9 r^ 1)/30 £ (0) - ^woo/? Ceo)- (54) 



'(0) 

Using then the first order Riemann tensor 

R^^p = ^ r ?i)At/3 ~~ ^r"^, (55) 

it reduces to 

+ <9or" )/30 C(o) = (^or")/3o - ^r^oo) £ (0) - (56) 



" MP 
dn 2 



'(0) 

Canceling doT^pQ on both sides, we get 



dn 2 

where 



(1) — a„r a „„ t@ 



Wi)OO^o), (57) 



i? 1)0 o = | ?f p (2 9o/i (1) po - dAuoo) • ( 58 ) 

Specializing now to the transverse-traceless coordinate system, where the components h^ p o 
vanish identically, we obtain 



0. (59) 



(0) 



Without loss of generality, we can take the solution to be = constant. In the linear 
approximation, therefore, in consonance with the fact that linear gravitational waves do not 
transport energy nor momentum, particles are not affected by linear gravitational wavesQ 

4.3 Second-Order Effects 

Up to second order, and already using the first order results, the geodesic deviation equation 
TED reads 



ds 2 



~ + r fi)70 r U)/?0 £<fi) + <9or" 2)/ 3o C( ) - ^(2)00y9 ^fo)- (60) 



'(0) 

Substituting the curvature tensor 



-^(2)00^ — do^wop dpFwoo + -^a)07^(i)0/3 ^a)^7^7i)00' (61) 



'For a detailed discussion of this point, see Ref. [10] 
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and considering that in transverse-traceless coordinates r7 1)00 = 0, we obtain 



d 2 e 

ds 2 



K 2 ) _ a pa cP 



dpT« oo e oy (62) 



'(0) 

Now, in transverse-traceless coordinates, the second order Christoffel connection is 

r (2)00 = 9ohf 2)0 - 7}d a h {2)00 . (63) 
The geodesic deviation equation reduces then to 

d 2 i%, (\ 



r/s 2 

U *(0) 



(2) - (^d a h (2)0 o ~ dpdoh^ ef 0) . (64) 



For definiteness, we consider a wave traveling in the z direction, in which case k p is 
given by Eq. (|41jh Let us then suppose two particles separated initially in the x direction 
by a distance £^ 0) , that is, 

eg) = (0,^,0,0). (65) 

Considering that, in the proper frame = ct, it is an easy task to verify that in this case 
the resulting equations of motion are 

^ 2 ^2) ^^2) p , fifi N 

at 2 5i 2 at 2 1 j 

The same result is obtained for two particles separated initially in the y direction. We 
consider now two particles separated initially in the z direction by a distance £ ( z 0) , that is, 

£f o) = (0,0,0,£f o) ). (67) 
In this case, the geodesic deviation equation (|64|) yields 

but 



1 <9 2 £ 2 

~^2~q]T~ = ( d odzh (2)z0 - \d z d z h {2)Q o) £ ( 2 0) . (69) 

This means that a gravitational wave does not produce movement orthogonal to the direc- 
tion of propagation. In other words, it is not an orthogonal, but a longitudinal wave. Notice 
that, in the second order, the two degrees of freedom are represented by h^zO = h^)0z and 

^(2)00- 

Considering that a detector on Earth will always be at large distances from the wave 
source, we use for h^uv the dominant solution (|45p . Furthermore, taking into account the 
arbitrariness of the wave vector K p , we can choose it in such a way that Kq = K\ = K<i = 0. 
In this case, the geodesic deviation equation ([69]) reduces to 

^L=^^. MM t-z/c)]. (70) 
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Although the wave amplitude decreases with distance, it can be assumed to be constant in 
the region of the experience. Accordingly, we write 



9 2 e 



2cz 

(2) 



1 



at 2 



£f 0) r (2) co 2 s'm[2(uit - z/X)], 



where 



(2) 



^ Z 



(71) 



(72) 



represents the wave amplitude at the region of the experience, with A = c/co the reduced 
wavelength. 

Observe that now the origin of the coordinate z is completely arbitrary. We can then 
choose one of the particles to be at z = 0, in which case z will represent the position of 
the second particle. Assuming that the particles are initially (t = 0) at rest, the solution is 
found to be 



£(0)^(2) 



sin[2(wt - z/X)} - 2ujtcos[2z/X] + sm[2z/A] 



(73) 



For gravitational waves with wavelength much larger than the particle separation (A S> z 
the solution becomes 

£(0)^(2) 



f'- 



sin(2a;i) — 2ut 



(74) 



(2) 16 

When a gravitational wave reaches two particles separated by a distance £^ 0) in the direction 
of the propagation, the distance between them will oscillate with frequency 2oj, and will 
grow linearly with time with a velocity 



6(0)^(2) w 



(75) 



This behavior is the result of tidal forces produced by the passage of a gravitational wave. 



5 Final Remarks 

Whenever use is made of a perturbation scheme, one forcibly ends up with a linear wave- 
equation. There is a widespread belief that gravitational waves can be approximately de- 
scribed by the solution of this linear wave equation. This assumption, however, is not 
justified. To understand it, let us make a comparison with gauge fields. As is well known, 
the gauge field of Chromodynamics must be nonlinear to transport color charge. Conversely, 
since electromagnetic waves are essentially linear, they are unable to transport their own 
source, that is, electric charge. Observe that, even though electromagnetic waves are lin- 
ear, they do transport energy and momentum. This is possible because neither energy 
nor momentum is source of the electromagnetic field. As such, the energy-momentum cur- 
rent does not enter the electromagnetic field equation, and consequently its linearity does 
not restrict the energy- momentum current to be linear. Differently from electromagnetic 
waves, however, in order to transport energy and momentum (the source of gravitation), a 
gravitational wave must necessarily be nonlinear. 
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If we accept that the first-order equations are fully correct up to that order, and not just 
nearly correct as is usually assumed in the mixed approach, we arrive at the inexorable result 
that linear gravitational waves transport neither energy nor momentum. As a consequence, 
they are unable to produce any effect on free particles. One may wonder why the first-order 
gravitational wave, which has a non-vanishing curvature tensor, produces no effects on free 
particles. To get some insight on this question, let us consider the following points. First, 
recall that the energy- momentum current is at least quadratic in the field variables H Since 
it appears explicitly in the gravitational field equations, in the linear approximation it is 
restrict to be linear, and consequently to vanish. This means that the energy-momentum 
density of any linear spacetime configuration must vanish. A non-vanishing energy density 
can only appear in orders higher than one. This does not mean that the first-order gravita- 
tional field is physically meaningless. In fact, at the second order it will appear multiplied 
by itself, giving rise to nonlinear field configurations with non- vanishing energy-momentum 
density. Second, observe that the components of the Riemann tensor are not physically 
meaningful in the sense that they are different in different coordinate systems. For exam- 
ple, starting with the "electric components" Riojo of the Riemann tensor, through a general 
coordinate transformation one can get non- vanishing "magnetic components" Riojk- By in- 
specting only the components, therefore, it is not possible to know whether they represent 
a true gravitomagnetic field, or just effects of coordinates. In order to get their physical 
meaning, one needs to inspect the invariants constructed out of the Riemann tensor compo- 
nents]^ Now, as a simple calculation shows, all invariants constructed out of the first-order 
Riemann tensor of the linear gravitational waves vanish identically [13] • This includes the 
scalar curvature, the Kretschmann invariant, and the pseudo-scalar invariant. Consider- 
ing that these invariants are proportional to the mass and angular momentum of the field 
configuration, we can conclude that the transverse components of the first-order wave are 
empty of physical meaning as no mass nor angular momentum (or helicity in the massless 
case) can be associated with them. 

Motivated by the above results, we have then obtained the second-order solution to 
the gravitational field equations, which might represent a physical gravitational wave. Its 
amplitude depends explicitly on the frequency of the wave — a property typically related 
to nonlinearity. In contrast to the linear wave, the second-order wave is able to transport 
energy and momentum. This fact becomes evident if we notice that, at second order, 
the source energy-momentum tensor is conserved only in the covariant sense; namely, it is 
not really conserved [15]. This means that, differently from what happens at first order, 
at second order a mechanical system can lose energy in the form of gravitational waves. 
Furthermore, although the first-order field is transverse and traceless, the second order is 
longitudinal. This property can be understood by remembering that gravitational waves 
are generated, and act on particles, through tidal effects, which arise from inhomogeneities 
in the gravitational field. The effects they produce on free particles are then found to 
consist of nonlinear oscillations along the direction of propagation. This is the signature a 
gravitational wave will leave in a detector, the effect to be looked for. 

A crucial point of these waves refer to their frequency. As is well known, the quadrupole 
radiation emitted by a source propagates with twice the frequency of the source (this has 

5 The classical energy of any wave, for example, is proportional to the squared amplitude. 
6 For a discussion of this point, see Ref. [T3], page 355. 



13 



to do with the tidal nature of the gravitational wave generation). However, the linear 
gravitational wave emerges from the perturbation scheme propagating with the same 
frequency of the source. To circumvent this problem, one has to artificially adjust by hands 
the wave frequency (see, for example, Ref. [2], page 105). On the other hand, owing to its 
nonlinear nature, the second order gravitational wave hf^) naturally emerges propagating 
with a frequency which is twice the source frequency, with the factor "2" coming from the 
fact that h 1 ^) depends quadratically on h?X. This is in agreement with the quadrupole 
radiation property, as well as with the tidal nature of the gravitational radiation, and is 
a clear indication that it is not the first-order, but the second-order wave that represents 
the quadrupole (physical) gravitational wave. It is also important to observe that, due 
to an explicit additional linear dependence on the source distance r, the amplitude of the 
dominant part of h^) is found to fall off as 1/r. Contrary to the usual belief, therefore, 
which presupposes that second-order effects fall off as 1/r 2 , second-order effects are not 
necessarily neglectful at large distances from the source. 

It is important to remark finally that, according to Birkhoff's theorem^ any spheri- 
cal source produces a time-independent gravitational field outside it. As a consequence, 
no spherically symmetric longitudinal gravitational waves can exist. However, due to the 
explicit dependence of the amplitude coefficient ([55]) on the wave number, we see that the 
nonlinear gravitational wave considered here will never be spherically symmetric. The usual 
restrictions imposed by Birkhoff's theorem on longitudinal gravitational waves, therefore, 
do not apply to the present case of longitudinal gravitational waves. 
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